Coupled Fluctuations near Critical Wetting 
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Recent work on the complete wetting transition has emphasized the role played by the coupling 
of fluctuations of the order parameter at the wall and at the depinning fluid interface. Extending 
this approach to the wetting transition itself we predict a novel crossover effect associated with 
the decoupling of fluctuations as the temperature is lowered towards the transition temperature 
Tw ■ Using this we are able to reanalyse recent Monte-Carlo simulation studies and extract a value 
ui(Tw) — 0.8 at -Sj- ~ 0.9 in very good agreement with long standing theoretical predictions. 
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A long standing controversy in the study of phase 
transitions at interfaces concerns the nature of the 
continuous wetting transition in three dimensional sys- 
tems with short ranged forces which corresponds to 
the marginal dimensionality 0~|^J^§1- Renormalization 
group (RG) analyses of simple capillary wave models 
predict [^J^] strong non-universality for critical exponents 
and amplitudes dependent on the value of the 'so called' 
wetting parameter lo(T) at the transition temperature 
Tw- However reliable estimates for lo(T) appropriate 
to the Ising model are significantly larger than values 
fitted to extensive Monte-Carlo simulation data [|| which 
only reveal small deviations from mean field (MF) 
theory [Q. Recently it has been suggested ^ that 
the transition in the Ising model is actually fluctuation 
induced (weakly) first-order hinting that the fitted values 
HQ are unreliable. However no quantitative analysis of 
the simulation data was given and the question, what 
value of w(T) is consistent with simulation studies of the 
Ising model remains unanswered. 

In this paper we extend our recent analysis of 
fluctuation effects at the complete wetting transition 
based on a 'two field' effective Hamiltonian ^[hjh] 
[|9|-|l"l| to the problem of the wetting transition discussed 
above. This approach predicts a novel crossover effect 
as the temperature is lowered towards Tw which is not 
predicted by simple capillary wave models (or their 
extensions incorporating a position dependent stiffness 
coefficient M) and is associated with the decoupling of 
fluctuations at the wall and at the fluid interface. This 
behaviour is in qualitative agreement with recent Monte- 
Carlo simulation studies of finite size (FS) effects in thin 
magnetic films [|l2| which have already been shown to be 
consistent with the two field theory deep in the complete 
wetting regime |j| [ll] . Using this prediction we are able 
to extract a value for the capillary parameter lu(T) at 
(or very close to) the transition temperature Tw whilst 
avoiding the issue of whether the transition is weakly 
first or second-order. We find cu{T w ) ^ 0.8 at ^ ~ 0.9 
which is in very good agreement with the most recent 
series expansion estimate at this temperature ||]. This 
is close to the universal critical value of u [Q which has 
been the long standing expectation for the magnitude of 
the capillary parameter appropriate to the Ising model 



To begin we make some remarks about the simulation 
studies used to extract this value of lu. Binder, Landau 



and Ferrenberg (BLF) |12| consider a thin Ising film 
(thickness D lattice spacings) with competing surface 
fields Hi = —Hrj. For this geometry there are 
a number of theoretical predictions for the way in 
which lcngthscales associated with wetting phenomena 
determine the nature of phase coexistence and criticality 
in the confined system Jl^,|lJ] . These are all confirmed by 
the BLF simulation studies. In particular BLF observe 
the predicted symmetry broken phase (in which wetting 
films are adsorbed at each wall) for temperatures T < 
T C (D) (with T C (D) <T W @) and a 'soft mode' phase 
Jl3| in the temperature window Tc > T > Tw (with Tc 
the bulk critical temperature) where an upspin-downspin 
interface sits on average at the centre of the thin film 
and whose fluctuations are controlled by an exponentially 
large correlation length 



£|| ~ exp(^) 



(1) 



Here k is the inverse bulk true correlation length and 
9 = 8(T,Hi) is a non-universal critical amplitude, the 
temperature dependence of which plays a crucial role in 
our analysis. Thus the qualitative features of both theory 
and simulation are in good agreement. 

The same is not true of the earlier simulations of 
Binder, Landau and co-workers |(| which studied wetting 
in a thin Ising film with equal surface fields H± = Hjj. 
This is the standard geometry modelling fluid adsorption 
between parallel plates where interest usually focuses on 
the phenomenon of capillary condensation corresponding 
to the coexistence of gas-like (downspin) and liquid-like 
(upspin) states at a FS shifted value of the bulk ordering 
field H. The observation of phases with wetting layers 
adsorbed at each wall for zero bulk field H = has 
been previously interpreted as corresponding to meta- 
stable states since the equilibrium must be a liquid-like 
capillary condensate. Such non-equilibrium states are 
observed above and below Tw- Indeed measurement of 
the wall susceptibility \i over this region forms the basis 
for locating Tw H and the method of extracting the 
value wet — 0.3 M which recall is much smaller than 
expected. This method of analysing the simulation data 
assumes that the behaviour of \i (and other response 



1 



functions) in the parallel plate geometry is similar to 
the thermodynamic semi-infinite limit. However this 
is not the case. We have considered the MF Landau 
type model for adsorption in the parallel plate geometry 
with identical surface fields and found that the 

wetting spinodal temperature Tyy , which marks the 
highest temperature (in zero bulk field) for which it 
is possible to have a meta-stable state with adsorbed 
wetting films, lies very close to and below Tw- Thus 
theory and simulation for this F S geometry are not even 
in qualitative agreement. It is therefore not surprising 
that the detailed quantative behaviour of the response 
functions do not yield the desired critical exponents. As 
mentioned above these problems do not arise for the 
FS geometry with reversed surface fields where the MF 
model and simulation results are in good qualitative 
agreement. For this case we can be fairly sure that 
the simulations are probing equilibrium behaviour which 
we hope to quantitatively understand using effective 
Hamiltonian theory. 

To continue we recall some pertinent ideas in the 
development of effective Hamiltonian models of wetting 
transitions. The standard capillary wave model H[l(y)] 
has the form |l],|| 



H[Ky)} 



dy 



£a/j(T) 



(yif+W(l(y)) 



(2) 



where T, a p(T) is the stiffness coefficient of the fluid 
interface (seperating bulk a and (5 phases) which unbinds 
from the wall and whose position is described by the 
collective coordinate l(y). For systems with short ranged 
forces the binding potential W(l) is taken to have the 
form |,| 



W(l) =hl + 



ac 



-kI 



be 



-2k/ 



(3) 



where k as mentioned above is the inverse bulk cor- 
relation length of the adsorbed {(3) phase |5]j|] and h 
is proportional to the bulk ordering field. At mean 
field level critical wetting occurs at h = a(T w lF ) = 
provided b > 0. Similarly the complete wetting transition 
occurs in the limit of vanishing bulk field h — > + for 
Tn > T > T$ F where a(T) > 0. RG theory based on 
(||) and ([}]) predicts critical exponents and amplitudes 
which are sensitive to the wetting parameter 



k B Tn 2 
47r£ a/3 



(4) 



(Ub being Boltzmann's constant). For values u> < 2 the 
phase boundary for critical wetting remains a(Tw) — 
i. e. Tw = T W IF but the critical exponents are very 
different from MF theory. For example along the isobar 
h = + the transverse correlation length diverges (as 
T — » T^) with an exponent v\\ = (\/2 — \/^)~ 2 f° r i < 
lu < 2 The implications for complete wetting are 

less dramatic — only critical amplitudes are sensitive to 
ui whilst exponents retain their MF values O. 



The discrepancy between these predictions and the 
older Monte-Carlo simulations HQ led Fisher and Jin 
Hpl| to reassess previous derivations of H[l(y)] conclud- 
ing that the stiffness coefficient should be replaced with 
a position dependent term 



E a/ 3(T) + ae- Kl - qle 



-2kI 



(5) 



although the binding potential W(l) is essentially correct. 
The important term in ([5]) is the next to leading order 
exponential decay which is negative (q > 0) at the MF 
phase boundary h = a{T^ F ) = 0. The coefficient a 
is proportional to a. When this position dependence 
is taken into account in RG calculations the wetting 
transition is driven first order for sufficiently small values 
of lo < lj* where the tricritical value to* is expected to be 
of order unity. Fisher and Jin estimate that the transition 
appropriate to the Ising model is very weakly first order 
(the correlation length is enormous at the transition) and 
that Tw is very close to T W * F ■ 

The concept of a position dependent stiffness co- 
efficient was forwarded independently by Parry and 
Evans ||I6| who pointed out that Hamiltonians of the 
form (0) could not describe next to leading order 
singularitcs of correlation functions at the complete 
wetting transition (which are known to exist from exact 
statistical mechanical sum rules). Unfortunately the 
position dependence of S(^) explicitly derived by Fisher 
and Jin using crossing and integral criteria is not of the 
type required by Parry and Evans fl6|| to satisfy full 
thermodynamic consistency. One way around this is 
to introduce a two field effective Hamiltonian H 2 [l\,l2\ 

M 

H 2 [W,l 2 ] = J dy[~SV(M2)VZ„. Vl v 

+U(h) + W {2) (l2-h)} (6) 

which models the coupling of fluctuations at the wall and 
af3 interface. The Hamiltonian (^) may be derived from 
an underlying 'microscopic' Landau-Ginzburg- Wilson 
functional Hlgw [ m ( r )] using a double crossing criteria in 
which h(y) and ^(y) are collective coordinates denoting 
surfaces of fixed magnetization and m 2 respectively. 
In the approach to the complete wetting transition the 
collective coordinate l 2 unbinds from the wall whilst l\ 
remains bound. The position dependence of the stiffness 
elements ^^(hih) provide a very elegant explanation 
of the correlation function singularites which single field 
Hamiltonians fail to describe. In particular it is the off 
diagonal elements 212(^21) ~ atd 2 \e~ td ' lr (with l 2 \ — 
h — h) which provide the dominant exponential decay 
essential for thermodynamic consistency. The term En 
is essentially position independent and may be identified 
with the stiffness of the wall-/3 interface T, w p whilst 
^22(^21) is very similar to the Fisher- Jin stiffness (||). The 
presence of coupled fluctuations has a rather profound 
effect on the critical behaviour at complete wetting. RG 
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calculations [p|]TT| predict that critical amplitudes are no 
longer determined by ui but by the renormalized quantity 



LU = LU + 



i + (Aie 



W0) 



(7) 



where tup 



4ttS„, 



and t; w p is the (finite) correlation 



length at the w/3 interface which may be related to 
T, w p and the curvature of U(l\). Ai is the momentum 
cut-off for the bound surface l\. The distinction 
between momentum cut-offs for the surfaces of fixed 
magnetization amd m£ was not addressed in our 
earlier discussion of coupling effects deep in the complete 
wetting regime. This will play an important role in our 
treatment of the crossover to critical wetting. 

Consider for example the effective Hamiltonian (||) 
with cut-off A. In standard interpretation the range 
of wavevectors allowed in the Fourier decomposition of 



i(y) is < Q < A where A <C y rfr corresponding 
to lengthscales much greater than the bulk correlation 
length [^7|. The same idea applies to the two field 
Hamiltonian. In fact it is easy to establish |lqj that 



the local cut-off Ai must satisfy Ai <C y otherwise 

this picture of fluctuations breaks down. This is a 
weak inequality and we expect Ai to vanish faster than 
^jYi w p. Of course it does not mean that fluctuations of 
the order parameter m(r) with wavevectors greater than 
Ai do not exist rather that only for sufficiently small 
wavevectors are they well described by an interfacial- 
like collective coordinate h{y) which couples to ^(y). 
Fluctuations in m(r) with wavevectors greater than Ai 
are not interfacial-like and have to be added to the model 
in some other way []l8|| . Coupling between these modes 
and I2 (y) does not lead to renormalization of the wetting 
parameter. 

We now apply these ideas to the wetting transition 
starting in the complete wetting regime and ask how the 
critical amplitudes describing the transition vary as the 
temperature is lowered towards T w tF ■ The important 
features of the Hamiltonian in this limit are 

(i) the cancellation of the leading order exponential 
decay in W%(l2i) and S22G21) similar to that indicated in 
(||) and (||) . Such behaviour could have been anticipated 
from the simpler capillary wave model and the Fisher- Jin 
theory. 

(ii) the vanishing of the local stiffness Y* w p ~ a 2 and 
hence the momentum cut-off Ai. These are new effects 
specific to the two field model. 

The important behaviour related to (ii) corresponds to 
a decoupling of fluctuations in the order parameter m(r) 
at the wall and a/3 interface. Within our theory all these 
features, (i) and (ii), are associated with the flattening 
of the magnetization profile near the wall as T — > T^ F . 
Whilst we were initially worried that this was an artifact 
of our model inspection of the simulation results for the 
magnetization p2| appears to confirm this very close to 
the observed T w . Exactly at T = T$ F the two field 



model is essentially identical to that of Fisher and Jin so 
repeating their argument § we predict that the wetting 
transition is fluctuation-induced (very weakly) first-order 
provided lu is not too big. However because the tricritical 
value lu* is not known very accurately the transition in 
the Ising model may still be second-order. Even if the 
transition was first-order it is very unlikely that this could 
be seen directly in simulations |ll[| . 

The new prediction of this analysis is a novel crossover 
effect associated with the decoupling of fluctuations. 
Deep in the complete wetting regime the coupling of 
fluctuations described by @ increases the effective value 
of the capillary parameter but this effect vanishes as the 
temperature is lowered to the wetting temperature (or 
more accurately T^ F ). This has important implications 
for the temperature dependence of critical amplitudes. 

In Fig. 1 we plot the values of the critical amplitude 
9 taken from the susceptibility measurements of BLF. 
At high temperatures T ^> Tw the observed value of 
the critical amplitude is bigger than the capillary wave 
prediction Ocw = 1 + % indicative that to should be 
replaced with To as predicted by the two field theory 
§,0. Recall the value of 9 at ^- ~ 0.96 is consistent 
with the estimate based on (^ which predicts an 
increment A9 ~ 0.3 JjHljJl. As the temperature is reduced 
there is clear evidence that 9 decreases consistent with 
the prediction discussed here. The extrapolated value 



lim 



(8) 



should be equal to I + ^ (if the transition is second order) 
and very close to this if it is (fluctuation induced) weakly 
first order. Using linear and cubic fits to extrapolate to 
the wetting temperature ^7 ~ 0.9 we obtain 8 + ~ 1.4 
which implies that lu(T) ~ 0.8 at this temperature. This 
is very close to the series expansion result of Fisher and 
Wen ||. Importantly there is no indication that the 'old' 
fitted values wg t — 0.3 and Ufu ~ || are appropriate. 
Such values of the capillary parameter would imply 9 + ~ 
1.15 or 9 + ~ 1 (corresponding to MF theory) which are 
totally inconsistent with the new BLF data. 
Four points that are worth emphasizing are 

(a) The method avoids the issue of whether the 
transition is first or second-order. Even if the transition 
is fluctuation induced first-order the wetting temperature 
Tw is very close to T^ F and coupling effects may be 
ignored when 9 is extrapolated to Tw- 

(b) We have carefully checked that the simulation and 
theory on the nature of FS effects in the parallel plate 
geometry are in qualitative agreement. We are confident 
that the simulations are probing equilibrium behaviour. 

(c) The data used to extract the critical amplitude 9 
are taken from the susceptibility measurements of BLF 
which are local to the fluctuating interface. Specifically 
they involve the susceptibilites Xn and x„„ centred at 
the middle of the thin film geometry where the interface 
lies on average. Older simulations H had relied on 
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FIG. 1. A cubic fit of the simulation results for 0(T, ■ ■ •). 
A linear fit to the lower four diamonds is also shown. The 
extrapolated values 6 + at -|tj ~ 0.9 are indicated by stars. 
Apart from the data point at ^ = 0.916 the error is within 
the symbol. As the temperature rises from TV the value of 6 
increases. 



measurements of the wall susceptibilites \i an d Xn 
which require assumption of scaling methods to analyse. 

(d) The correlation length £|| in the soft mode phase is 
very large and easily satisfies the Ginzburg criteria p9[ . 

In summary we have used a two field Hamiltonian 
model to predict a novel decoupling effect at the three 
dimensional wetting transition which has implications 
for the temperature dependence of critical amplitudes. 
On the basis of this we have reanalyzed recent Monte- 
Carlo simulations and extracted a value of the wetting 
parameter which is very close to long standing theoretical 
predictions. 
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